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With the Finite temperature Density Matrix Renormalization Group method (FT- 
DMRG), we depeloped a method to calculate thermo-dynamical quantities and the con- 
ductance of a quantum dot system. Conductance is written by the local density of states on 
the dot. The density of states is calculated with the numerical analytic continuation from the 
thermal Green's function which is obtained directly from the FT-DMRG. Typical Kondo be- 
haviors in the quantum dot system are observed conveniently by comparing the conductance 
with the magnetic and charge susceptibilities: Coulomb oscillation peaks and the unitarity 
limit. We discuss advantage of this method compared with others. 
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1. Introduction 

Recently, the Kondo effect is an attractive topic in quantum dots (QD). At high tem- 
peratures, the conductance has peaks as a function of the gate voltage when the number of 
electrons in the dot changes, i.e. coulomb oscillation peaks. As a result of the Kondo effect, the 
conductance becomes larger up to the unitarity limit (2e 2 /h) when the number of electrons 
is odd. 1 Since a dot with an odd number of electrons plays a role of a magnetic impurity, the 
conduction electrons in the source and drain leads tend to screen the local moment by the 
Kondo effect. What is different from magnetic impurities is that the number of electrons is 
tunable from a magnetic region to a non-magnetic region by changing the gate voltage. Ad- 
ditionally, in a lateral QD the strength of the tunnel coupling between the dot and the leads 
also changes by changing the gate voltage. In this sense, not only temperature dependence but 
also dependences of impurity parameters are important to understand the quantum transport 
through QDs. 

For theoretical study of the conductance through QDs in a wide range of parameters, the 
Finite Temperature Density Matrix Renormalization Group Method (FT-DMRG) 2 ' 3 is used 
in this paper. FT-DMRG has the advantage that we can re-use the result for the bulk part for 
different impurities or dots. Strictly speaking, what we should do for various impurities is to 
average a quantum transfer matrix of the impurity site with the left and right eigen vectors 
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which are obtained in the bulk-part calculations. Additional advantage of the FT-DMRG 
is that not only static properties but also dynamical properties of the impurity site can be 
calculated. However, FT-DMRG has not been applied to quantum dot systems, to the best 
of our knowledge. One of the reasons may be the difficulty of the analytic continuation which 
will be discussed in the next paragraph. For the application of the FT-DMRG and White's 
DMRG method, 4 the system must be one-dimensional, but this is not a serious problem for 
QD because we can use the model of QD as a impurity between semi-infinite tight-binding 
chains. 

On the other hand, Wilson's Numerical Renormalization Group (NRG) 5-7 method is one 
of the most powerful methods to deal with the Kondo effect, because the NRG is designed for 
lower temperature than the Kondo temperature which is exponentially small by focusing on 
the low energy excitations in logarithmic scale. By the FT-DMRG it is difficult to reach that 
small temperature, while the high energy excitations are treated correctly. This limitation of 
the lowest temperature originates from the numerical errors of the algorithm. The numerical 
analytic continuation which is used to obtain Green's functions in frequency space from ther- 
mal Green's functions is the main source of the difficulty to decrease temperatures. To avoid 
the difficulty we use relatively large tunnel coupling and coulomb energy in order to make the 
Kondo temperature high. Fortunately the conductance which can be obtained from Green's 
function after the numerical analytic continuation has a better accuracy than the Green's 
function itself. 

In the following sections, we will show that in a certain parameter range conductance is 
obtained reliably from the FT-DMRG and show the typical Kondo effect in QD: the Coulomb 
oscillation peaks and the unitarity limit. In §3 we will discuss the parameter dependence of 
the Kondo temperatures in the symmetric case where the Kondo temperature has the lowest 
value as a function of the gate voltage. In §4 we will show that there is the set of parameters 
which meets the condition and plot the typical Kondo effect in QD from the point of view 
not only of conductance but also of thermodynamic quantities. In the next §2, we start from 
definition of the model Hamiltonian and demonstration of how to calculate conductance by 
looking at the non-interacting case. 

2. Hamiltonian and the conductance in the non-interacting case 

FT-DMRG can be applied to various impurities not restricted to a single dot. However 
in this paper we will concentrate on the Kondo effect in the single dot. For this purpose, we 
will restrict ourselves to the simplest Hamiltonian, which is a tight-binding chain with a dot 
which has a single local orbital, 
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where the site i = indicates the dot which has on-site Coulomb energy U and e^ corresponds 
to the gate voltage. This Hamiltonian is sometimes called the Anderson- Wolff type impurity 8 
(Fig.l). Magnetic field h and Chemical potential /j, are applied on the whole system. We use 
t as the units of energies. This simplest model has been used by many authors for the study 
of QD, because it includes the essential physics of the Kondo effect observed in experiments. 

With this Hamiltonian, we will calculate thermodynamic quantities and a conductance 
in h = n = case. Concerning the former, magnetic susceptibility (ximp = d(s^ =0 ) / dh\h=o) 
and charge susceptibility (ximp,c = <9(nj = o)/<9/i| M =o) are calculated to study behaviors of spin 
and charge fluctuations of the impurity site. By this definition the units for the magnetic 
susceptibility are (g/iB) 2 - It should be emphasized that basic properties of these quantities are 
the same for the s-wave impurity- Anderson model. This is easy to show, because this simplest 
Hamiltonian for QD can be mapped to two semi-infinite chains by using the liner combination 
such as -j=(a a ± Qo-); even-parity and odd-parity channels. This situation is completely the 
same with the s-wave impurity-Anderson model attached to one-dimensional chain. On the 
other hand, transport properties of the two models are very different. However, in the present 
paper we will restrict ourselves to the model given by eq.( 1). 

The latter, a conductance, can be calculated as follows. Thermal Green's function in 
imaginary time is obtained from the FT-DMRG by calculating an expected value with respect 
to the eigen vectors having the largest eigen vectors. In order to calculate the Green's function 
in frequency space from the thermal Green's function, we used the Maximum entropy method 
(MEM) 9 and the Pade approximation 10 as numerical analytic continuation. Since the MEM 
is basically a method to deal with statistical data, this method has been used with QMC 
to obtain the DOS of the QD 12 . However, the MEM was also applied to the FT-DMRG 
and it has been shown that the MEM is more reliable than the Pade approximation. 11 After 
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Fig. 2. Conductance as a function of the temperature for U — at various v calculated by the FT- 
DMRG with the finite algorithm. The data are obtained with the MEM and the Pade from Go(r). 
Trotter number M < 200 and the residual bases m < 145 are used. The solid lines are analytic 
results: p(e) = 



■nv 2 ' 4— v~ 2 (2 — v~ 2 )e 2 ' 



the numerical analytic continuation, we can calculate the conductance from the following 
equation, 12 ' 13 when the hopping terms from the impurity site to the left and right are the 
same. 
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where Go( e ) is the Green's function of the impurity site i = and A L / R is the local DOS of 
the next sites to the impurity site; Al = Ar = A = v 2 ^jAt 2 — e 2 /(2£ 2 ) and / is the Fermi 
distribution function. It should be noted that A has some e dependence in order to calculate 
the conductance at high temperatures which are comparable to the band width. It should 
be also noted that at zero temperature the conductance follows a simple expression 14 by the 
Freedel sum rule: 
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where (rid) is the number of electrons at zero temperature and is defined as the difference 
between the total number of electrons iVtot and the system size iV i.e. (rid) = (-^tot) — {N — 1). 
Since (rid) is equal to 1 in the symmetric case (e^ = —{7/2), we obtain quantized conductance 
g = 2e 2 /h at zero temperature. 

Before starting discussions on the Kondo effect of QD, we first demonstrate reliability of 
the conductance obtained by the FT-DMRG with the numerical analytic continuation for the 
non- interacting case (e^ = —U/2 = 0). In Fig. 2, we compare the results obtained by the 
Pade and the MEM with the exact ones (the solid lines in Fig. 2). Conductances calculated 
with using these two methods show the tendency that errors become more significant closer to 
the unitarity limit: g = 2e 2 /h. In fact, at lowest temperatures, some data are larger than the 
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Fig. 3. The schematic diagram in the symmetric Anderson model, reproduced from Fig. 12 in 6 

unitarity limit and these unphysical data are out of the range of the panel and not displayed. 
In conclusion we may say that the numerical results of the conductance are reliable when the 
conductance is away from the unitarity limit, except for the v = 1 case where the local DOS 
at the impurity site becomes singular at the band edge. 

3. Kondo singlet and Unitarity limit for the symmetric case 

Since FT-DMRG is more reliable at higher temperatures, we should arrange parameters 
in this Hamiltonian so that the Kondo temperature is sufficiently high. Since the Kondo 
temperature as a function of the gate voltage has the lowest value in the symmetric case 
(e^ = —U/2) , the symmetric case is the most difficult case to see the Kondo effect. In this 
section, therefore, we restrict ourselves to the symmetric case and discuss the dependence on 
the other parameters, U and v. In the symmetric case, the impurity susceptibility shows a 
Curie low (T% ~ 1/4 for 1/2 spin) in the temperature range where the local moment is well 
defined in the sense that the coupling to neighboring sites is negligible. However, due to the 
Kondo effect the local moment is screened (x = Const.) at the low temperature limit: the 
strong-coupling regime in terms of the scaling theory. In contrast to the Kondo model or the s- 
d model, the Anderson model has another high temperature limit (T\ ~ 1/8 at T> irv 2 , U): 
the free-orbital regime . 6 The authors of 6 wrote a schematic diagram like Fig. 3. In this 
schematic diagram we find that there is a direct cross-over from the free-orbital regime to the 
strong-coupling regime in the region irv 2 > U. This strong-coupling regime is continuously 
connected to the non interacting limit U = 0. In fact, magnetic susceptibilities Tximp for a 
typical case of U = 2 (Fig. 4) obtained from the FT-DMRG show a good consistency with 
Fig. 3. Each lines start from the free-orbital regime at high temperatures. Especially v = 1 line 
shows a direct cross-over to the strong-coupling regime. On the other hand, v = 0.2 line shows 
clearly the local-moment regime, but the strong-coupling regime cannot be seen because the 
Kondo temperature Tk is lower than 0.02. 

Similarly, conductance also shows a consistency with the schematic diagram. The low 
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Fig. 4. Susceptibilities T%i mp /(c;/iB) 2 as the function of temperature for U = 2 and various v cal- 
culated by the FT-DMRG with the infinite algorithm. The truncation error is 1CP 3 ~ 10~ 4 , the 
maximum number of the residual bases m is less than 145 and the Trotter number M is less than 
200. 

(and high) temperature limit can be understood from eq. 3 (eq. 4) and the conductance g is 
equal to 2e 2 jh at zero temperature. In fact, conductance for a typical case U = 2 (Fig. 5) 
shows that each line starts from the free-orbital regime g ~ at high temperatures. Especially 
v = 1 line shows a direct cross-over to the strong-coupling regime g ~ 2e 2 //i. On the other 
hand, v = 0.2 line does not show clearly the local-moment regime in contrast to the magnetic 
susceptibility. The local-moment regime in conductance seems to be difficult to distinguish 
from the free-orbital regime if we identify the local-moment regime simply by vanishingly 
small g. This naive identification is consistent with the v = case where g = at any 
temperature. However v = 0.2 line in Fig. 5 shows a slight peak at T ~ 1, which may be 
on indication of the crossover between the free-orbital regime and the local-moment regime. 
Strictly speaking, this broad hump comes from the double-peak structure of the local DOS 
at w = e, e + U through the integration (eq. 3). Then, we find that the peak position of 
this conductance is around T ~ U/2 and it's height becomes zero at v = 0. It should be 
noted that the results based on Pade approximants show sometimes instabilities in Fig. 5 
(for example, v = 0.4 and T = 0.15), because numerical errors are accidentally accumulated 
in the iteration process of the Thiele's reciprocal method. 10 Because of the numerical analytic 
continuation, the conductance is more difficult quantity than the magnetic susceptibility. 
However, conductance obtained after integration in frequency space is fairly reliable except 
for some special cases as is seen in Fig. 5 

From Fig. 4 and 5 one can see that it is difficult to reach the strong coupling regime 
by the present method within the lowest temperatures when v is small. To determine the 
region where the strong coupling regime is reachable, we fix the temperature (T = 0.1) 
and plot the magnetic susceptibility (Fig. 6) and the conductance (Fig. 7) as a function 
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Fig. 5. Conductance as a function of the temperature for U — 2 with various v calculated by the 
FT-DMRG. The other conditions are the same as in figure 2. 







TX m atT=0.1 

contour — +- 




Fig. 6. The magnetic susceptibility Tximp as a function of U and v in the symmetric case at fixed 
T = 0.1. Three regimes are identified :the strong-coupling regime as Tximp ~ , the local- 
moment regime as Txi mp ~ 1/4 and the free-orbital regime around the origin, which is the region 
ttv 2 ,U <T. 

of v and U. Note that the z-axis of Fig. 6 increases downward. These two graphs show 
a similar behavior concerning the crossover from the local-moment regime to the strong- 
coupling regime. The difference between the two figures is concerned with the free-orbital 
regime (around irv 2 < T and U < T), because the local- moment regime and the free-orbital 
regime are hard to distinguish for the conductance. With the remark in mind, one can conclude 
that the cross-over behaviors from the local-moment regime to the strong-coupling regime are 
consistent between the two quantities. 

Finally, we try to estimate the cross-over temperature which is nothing but the Kondo 
temperature. At zero temperature, all parameter region goes to the strong-coupling regime 
(Tximp = and g = 2e 2 /h) under the symmetric condition. Physical aspects of the 
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Fig. 7. The conductance g as a function of U and v in the symmetric case at fixed T = 0.1. The units 
of g are 2e 2 /h. 

strong-coupling regime appear in the value of Ximp itself at zero temperature and in the 
way how g approaches to the unitarity limit. In U — ► oo limit, the susceptibility at 
zero temperature is identified as Ximp = l/(4Tk), where Tk means the Kondo tempera- 
ture. This Kondo temperature also determines the temperature dependence of conductance: 
g = (2e 2 /h) (1 — fg(— ) 2 + 0(T 4 )j. The latter equation comes from the temperature de- 
pendence of the Kondo resonant peak in the DOS. Figure 8 shows the Kondo temperature 
obtained from Ximp an d also from g (Inset) by using the definitions of the large U limit. Since 
numerical data of g do not have enough accuracy near the unitarity limit, it is hard to esti- 
mate the coefficient of T 2 in g. Therefore, we determine Tk simply from the temperature T 
where g = 0.75 with that equation. Tk thus estimated from the conductance is only a rough 
estimation. However, the two figures are consistent with each other. Strictly speaking, these 
formulas which are used to define Tk are exact in the U — ► oo limit where the Anderson model 
is mapped to the s-d model. We consider that Tk shown in Fig. 8 qualitatively expresses the 
cross-over temperature to the strong-coupling regime. In this sense, we may consider that 
qualitative figure of Fig. 3 are reproduced numerically by Fig. 8. 

As a conclusion of this section, we have observed behaviors in the strong-coupling regime 
both in magnetic susceptibility and conductance. Concerning the latter, there is a parame- 
ter range where the conductance reaches the unitarity limit within the lowest temperature 
determined by the stability of the numerical analytic continuation. It is concluded that the 
strong-coupling regime is easier to see when v 2 /U is large. However, to reproduce the typical 
Kondo behaviors which occur in the local-moment regime, another limitation (v 2 /U is small) 
is needed. In the next section, we will show that there is a region of parameters which sat- 
isfy the both limitations. By choosing a parameter in the region, one can observe that the 
conductance changes from the Coulomb oscillation regime to the unitarity limit. 
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Fig. 8. The Kondo temperature 7k estimated from the susceptibility Ximp at the lowest temperature 
T = 0.01. Inset: The Kondo temperature defined as 7k | g =o.75= \T | s= o.75, where T | 3= o.75 is 
defined as the temperature where the conductance g becomes 0.75(2e 2 //i). This estimated Kondo 
temperature is plotted only at high temperature region, because data of the conductance are 
available only T > 0.1 . 
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Fig. 9. Conductances as a function of the gate voltage V g for v = 1/V2, U = 4 at various temperatures 
calculated by the FT-DMRG with the MEM and the Pade . The other conditions for the numerical 
calculation are the same as in figure 2. Since td = —U/2+V g , V g = means the symmetric condition 
and V g = 2 means €d = 0. 

4. Coulomb oscillation for the asymmetric case 

In QD, the number of electrons is tunable by changing the gate voltage which corresponds 
to changing Vg in the asymmetric case (e<f = —U/2 + Vg). We stress again that FT-DMRG 
has the advantage to deal with the cases of different gate voltages at once. As an important 
character of the asymmetric case, Coulomb oscillation peaks appear in the conductance as 
a function of the gate voltage, because the conductance becomes non-zero only when the 
rid = and n^ = 1 configurations, or the n<\= 1 and n^ = 2 configurations, are degenerated 
, so-called the valence-fluctuation regime. On the other hand, it is easy to understand the 
behavior of a conductance at zero temperature, because the conductance is determined by 
rid (eq. 4). Especially in the case of {rid) ~ 1 5 the conductance reaches the unitarity limit. 
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As we mentioned in the previous section, there is numerical difficulty to see the Coulomb 
oscillation and the unitarity limit with the same U. Since the Coulomb oscillation peak comes 
from the double peaks of local DOS at u> = e^, e^ + U, the condition U 3> ttv 2 is required. We 
see clear coulomb oscillation peaks by putting v 2 /U smaller, but it makes difficult to see the 
strong-coupling regime within T > 0.1. Then, to compromise the conflicting requirements, we 
take U = 4 and v = l/v2 ~ 0.7, where U is slightly larger than ttv 2 = tt/2. The conductance 
for this set of parameters is shown in Fig. 9. In fact, Fig. 9 shows that the coulomb oscillation 
peaks at V g = ±2 are overlapping. However, this figure shows the typical Kondo effect in the 
QD: Coulomb oscillation peaks gradually move towards V g = at low temperatures and fill 
the valley up to the unitarity limit. 

Again, it is instructive to compare conductance and magnetic susceptibility. To make the 
situation simple, let us start with the atomic limit (v = 0). For the magnetic susceptibility in 
the atomic limit, there are several fixed points: 6 the free-orbital fixed point (Txi mp = 1/8), the 
valence-fluctuation fixed point (T^mp = 1/6), the local-moment fixed point (Tximp = 1/4) 
and the frozen- impurity fixed point (Tx; mp = 0). It is easy to calculate the susceptibility in 
the atomic limit: Txi mp =e t < 2 I 1 + 2e ^ + e t \\ . Actually, we can see the above 
mentioned regimes in the plot of this function for U = 4 as shown in Fig. 10(a) . Figure 
10(b) is Tximp obtained from the FT-DMRG for U = 4 and v = —?=. One can see that a 
developed moment in the local-moment regime, |V^| < U/2 and T < U, is very suppressed by 
the effect of v. The valence fluctuation regime, which is characteristic to the asymmetric case 
and is expected to correspond to the Coulomb oscillation, is not clearly seen due to relatively 
large value of v. The valence fluctuation regime may be characterized by enhanced charge 
fluctuations. Therefore, we plot the charge susceptibility in Fig. 11(b) . 

To understand the results we plot in Fig. 11(a) the charge susceptibility Tximp,c m the 
atomic limit v = for U = 4. In the v = case, it is straightforward to obtain: Tx\mp,c = 

t d + U ( fd 2e d + U \ / c d + U 2e d +U \ -2 

2e t M -)- 2e r + e t 1 + 2e t +e t \ . The value of Txi m p,c is 1/2 in the free- 



orbital regime and 1/4 in the valence fluctuation regime and zero in the other regimes. Then 
we find that sharply defined valence fluctuation regime at \V g \ ~ U/2, Fig. 11(a) , is suppressed 
by the finite v ,Fig. 11(b) . This behavior is the same as the magnetic susceptibility. However, 
there remain small structures which reflect the suppressed valence fluctuation regime and are 
thought to correspond to the Coulomb oscillation peaks. 

Finally, we discuss low temperature properties of QD by comparing Ximp (Fig- 12), Ximp,c 
(Fig. 13) and g (Fig. 9). When |V^| > 2, temperature dependences of these quantities show 
convergence already at T = 0.1. On the other hand, when the number of electrons is odd 
(|V^| < 2), the fixed point at zero temperature is far from the lowest temperature T = 0.1 
in the present calculations, as can be seen from the conductance which does not reach the 
unitarity limit at least in the symmetric case. This behavior is consistent with the fact that 
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Fig. 10. (a) Three dimensional plot of a spin susceptibility Tx; mp as a function of T and V g for 
U = 4 in the atomic limit v = 0. T%i mp = ° + ^ — 2 td +u ■ (b) Three dimensional plot 

2(l+2c T +c T ) 

of a spin susceptibility Tximp as a function of T and V^ for v = l/v2- Data are plotted at 
T = 0.1, 0.15, 0.3, 0.6, 1.2, 3, 6, 12, 30. 



Ximp has a maximum at the symmetric point (Fig. 12), because the inverse of Ximp(T = 0) 
is proportional to the cross-over temperature when |Vg| < 2, i.e. a large Ximp means a small 
cross-over temperature. On the other hand, the charge susceptibility Ximp,c(T = 0) and the 
conductance g(T = 0) reflect the number of electrons as the function of e^. Their behaviors 
are quite different when the number of electrons is odd. However, their basic behavior is easy 
to understand because the number of electrons changes like a step function with increasing 
ed- The reason why Ximp,c is not zero at V^ = in Fig. 13 is that v 2 /U is finite and not small. 
When v 2 /U is small, n^ becomes flatter around the symmetric point and Ximp,c becomes 
smaller. 

It can be concluded that we can observe the valence-fluctuation regime as visible peaks 
in the conductance and in the charge susceptibility, because the valence-fluctuation regime 
is sensitive to the conductance and also to the charge susceptibility when V g is swept at a 
fixed temperature. These peaks give rise to Coulomb oscillation peaks at higher temperatures. 
When the strong-coupling regime appears at lower temperatures, the charge susceptibility and 
the magnetic susceptibility are saturated, i.e. Ximp,c an d Ximp takes finite values. Since the 
cross-over temperature from the valence-fluctuation to the strong-coupling regime is thought 
to be higher than that at the symmetric point V g = 0, Coulomb oscillation peaks in the 



11/15 



J. Phys. Soc. Jpn. 



Full Paper 




0.5 
0.45 

0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 

0.05 





Fig. 11. (a) Three dimensional plot of a charge susceptibility Tx; mP]C as a function of T and V g for 
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Fig. 12. A spin susceptibility Ximp as a function of V g at various T for U = 4 and f = l/v2- 

conductance gradually move to V^ = with decreasing temperature as we see in Fig. 9. At 
last, the conductance at V g = should reach the unitarity limit, because the conductance 
at zero temperature follows eq. 4. However, at low temperatures where the conductance is 
close to the unitarity limit, there is the difficulty of the numerical analytic continuation as we 
mentioned in the previous section. 
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Fig. 13. A charge susceptibility Ximp,c as a function of V g with various T for U = 4 and v = l/v2- 

5. Summary 

We have developed a method to calculate thermo-dynamical quantities and the conduc- 
tance at various parameters simultaneously by using the FT-DMRG, because FT-DMRG can 
use the same result for the bulk part in the calculations of various v, U and e^. This is the 
advantage of the FT-DMRG; for example, the bulk-part calculation needs a large memory 
mainly in a Lanczos routine for the quantum transfer matrix, but the impurity-part calcu- 
lation does not need such a process. What we should do next is to calculate the average of 
physical quantities. Actually we have seen that cross-over behaviors between different regimes 
are observable not only in the conductance but also in the magnetic and charge susceptibilities 
by changing impurity parameters. Especially we can clearly see nice correspondence between 
the Unitarity limit and the strong-coupling regime and also between the Coulomb oscillation 
peaks and the valence-fluctuation regime. We may conclude that FT-DMRG is useful method 
to study the QD. 

Next we will discuss some technical aspect of the FT-DMRG. Since FT-DMRG is an 
approach where the numerical error accumulates as temperature is lowered, we can not use 
realistic parameters i.e. v and U are much smaller than the band width t, because in the 
parameter region the Kondo temperature becomes extremely small. In principle, we can sepa- 
rate up and down spin electrons in the bulk part calculation. By doing this, probably one can 
go to lower temperatures than the present result. The reason we did not use the algorithm, 
is that we will use the present program to more general system such as the case that the 
Coulomb interaction exists in the leads. The effect of many body interaction in the leads will 
be reported in a separate publication. 

In order to obtain reliable results at low temperatures with this method, there is another 
problem: the numerical analytic continuation. For example, as we mentioned already raw data 
of conductance sometimes become larger than the unitarity limit at T ^C Tk- The numerical 
error comes from two reasons; one is an error of the FT-DMRG in a thermal Green's function 
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Gq(t), and the other is an additional error coming from the numerical analytic continuation 
from Gq(t) to Go(e). The numerical errors of Gq(t) at r = 0,(3 correspond to truncation 
errors, because Gq(t = 0) is the number of electrons at the dot. The truncation errors are 
typically 10~ 3 ~ 10~ 4 in this calculation. Since G{t ~ /3/2) becomes smaller with decreasing 
temperature, the relative error becomes larger at lower temperatures. G(t ~ /3/2) has a 
stronger dependence both on the Trotter number M and on the maximum number of residual 
bases m. This is the difficulty concerning the numerical calculation of the thermal green 
function. In addition to this, the numerical analytic continuation generate additional errors, 
which appear even if Gq(t) calculated with a high accuracy from the explicit DOS is used. 15 
However, as we mentioned in §2 concerning the noninteracting case, this error becomes small 
if the conductance is away from the unitarity limit. It means that one can expect a better 
results for the conductance than the DOS itself, because even if the DOS has some error as 
the function of uj, it's integrated value like the conductance may show a better behavior. 
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